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' The shortest path between two truths in the real domain passes through the 

' complex domain. — Jacques Hadamard, The Mathematical Intelligencer 13 (1991) 

(N : 

^ I Consider the set of all Hamiltonians whose largest and smallest energy eigenvalues, 

' -E'max and -Emiru differ by a fixed energy oj. Given two quantum states, an initial state 

• lipi) and a final state \iPf), there exist many Hamiltonians H belonging to this set 
. under which evolves in time into IV'f)- Which Hamiltonian transforms the initial 

,— i, state to the final state in the least possible time r? For Hermitian Hamiltonians, r 

1^ I has a nonzero lower bound. However, among complex non-Hermitian PT-symmetric 

^ ' Hamiltonians satisfying the same energy constraint, r can be made arbitrarily small 

Ch ' without violating the time-energy uncertainty principle. The minimum value of r 

^ . can be made arbitrarily small because for PT-symmetric Hamiltonians the evolution 

O^! path from the vector \ipj) to the vector \^pF), as measured using the Hilbert-space 

I metric appropriate for this theory, can be made arbitrarily short. The mechanism 

' described here resembles the effect in general relativity in which two space-time 

• points can be made arbitrarily close if they are connected by a wormhole. This 

CN ■ result may have applications in quantum computing. 

00 ■ 
^ ■ 

25 ' I- INTRODUCTION 

oo , 

O ■ Interest in optimal time evolution dates back to the end of the seventeenth century, 

>■ , when the famous brachistochrone problem was solved almost simultaneously by Newton, 
I Leibniz, I'Hopital, and Jacob and Johann Bernoulli. The word brachistochrone is derived 
^ ' from Greek and means shortest time (of flight). The classical brachistochrone problem is 
stated as follows: A bead slides down a frictionless wire from point A to point 5 in a 
homogeneous gravitational field. What is the shape of the wire that minimises the time of 
flight of the bead? The solution to this problem is that the optimal (fastest) time evolution 
is achieved when the wire takes the shape of a cycloid, which is the curve that is traced out 
by a point on a wheel that is rolling on flat ground. 

In the past few years there has been much interest in the quantum brachistochrone prob- 
lem, which is formulated in a somewhat similar fashion: Consider two fixed quantum states, 
an initial state {ipi) and a final state {ipp) in a Hilbert space. We then consider the set of 
all Hamiltonians satisfying the energy constraint that the difference between the largest and 
smallest eigenvalues is a fixed energy u: -Emax — E^in = ^- Some of the Hamiltonians in this 
set allow the initial state {ipi) to evolve into the final state {ipp) in time t: 

|^^)=e-^*/«|^,). (1) 



2 



The quantum brachistochrone problem is to find the optimal Hamiltonian; that is, the 
Hamiltonian that accomphshes this evolution in the shortest possible time, which we denote 
by T. 

In this chapter we show that for Hermitian Hamiltonians the shortest evolution time r 
is a nonzero quantity whose size depends on the Hilbert-space distance between the fixed 
initial and final state vectors. However, for complex non-Hermitian Hamiltonians, the value 
of T can be made arbitrarily small. Thus, non-Hermitian Hamiltonians permit arbitrarily 
fast time evolution. 

Of course, a non-Hermitian Hamiltonian may be physically unrealistic because it may 
possess complex eigenvalues and/or it may generate nonunitary time evolution; that is, 
time evolution in which probability is not conserved. However, there is a special class 
of non-Hermitian Hamiltonians that are PT symmetric; that is, Hamiltonians that are 
invariant under combined space and time reflection. Although such Hamiltonians are not 
Hermitian in the Dirac sense, they do have entirely real spectra and give rise to unitary time 
evolution. Thus, such Hamiltonians define consistent and acceptable theories of quantum 
mechanics. We show in this chapter that if we use Hamiltonians of this type to solve the 
quantum brachistochrone problem, we can achieve arbitrarily fast time evolution without 
violating any principles of quantum mechanics. Thus, if it were possible to implement faster- 
than-Hermitian time evolution, then non-Hermitian Hamiltonians might have important 
applications in quantum computing. 

This chapter is organised as follows: In Sec. [ITl we introduce and describe PT quantum 
mechanics and explain how a Hamiltonian that is not Dirac Hermitian can still define a 
consistent theory of quantum mechanics. Then in Sec. IIIII we explain why complex clas- 
sical mechanics allows for faster-than-conventional time evolution. In Sec. [IV] we discuss 
the quantum brachistochrone for Hermitian Hamiltonians. Then, in Sec. |V] we extend the 
discussion in Sec. IIVI to Hamiltonians that are not Dirac Hermitian. In Sec. IVII we explain 
how it might be possible for a complex Hamiltonian to achieve faster-than-Hermitian time 
evolution. 

II. PT QUANTUM MECHANICS 

Based on the training that one receives in a traditional quantum mechanics course, one 
would expect a theory defined by a non-Hermitian Hamiltonian to be physically unacceptable 
because the energy levels would most likely be complex and the time evolution would most 
likely be nonunitary (not probability-conserving). However, theories defined by a special 
class of non-Hermitian Hamiltonians called PT-symmetric Hamiltonians can have positive 
real energy levels and can exhibit unitary time evolution. Such theories are consistent quan- 
tum theories. It may be possible to distinguish these theories experimentally from theories 
defined by Hermitian Hamiltonians because, in principle, non-Hermitian Hamiltonians can 
be used to generate arbitrarily fast time evolution. 

We use the following terminology in this chapter: By Hermitian, we mean Dirac Hermi- 
tian, where the Dirac Hermitian adjoint symbol ^ represents combined matrix transposition 
and complex conjugation. The parity operator P performs spatial reflection and thus in 
quantum mechanics it changes the sign of the position operator x and the momentum oper- 
ator p: PxP = —X and PpP = —p. Because the parity operator P is a reflection operator, 
its square is the unit operator: P^ = 1. The time-reversal operator T performs the time 
reflection t — —t, and thus it changes the sign of the momentum operator p, TpT = —p, 
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but it leaves the position operator invariant: TxT = x. The square of T is the unit operator 
= 1. We require that the operators P and T individually leave invariant the fundamen- 
tal Heisenberg algebra of quantum mechanics [x,p] = i. Thus, while P is a linear operator, 
we see that T must perform complex conjugation TzT = z*, and hence T is an antilinear 
operator. 



The first class of PT-symmetric quantum-mechanical Hamiltonians was introduced in 
1998 [gI]. Since then there have been many papers on this subject by a wide range of 
authors. There have also been three recent review articles 0, H, 29]. In Ref. it was 
discovered that even if a Hamiltonian is not Hermitian, its energy levels can be all real and 
positive so long as the eigenfunctions are symmetric under PT reflection. 

These new kinds of Hamiltonians are obtained by deforming ordinary Hermitian Hamil- 
tonians into the complex domain. The original class of PT-symmetric Hamiltonians that 
was proposed in Ref. has the form 

H = p^ + x^{ixy (e>0), (2) 

where e is a real deformation parameter. Two particularly interesting special cases are 
obtained by setting e = 1 to get H = + ix^ and by setting e = 2 to get H = p^ — x^. 
Surprisingly, these Hamiltonians have real, positive, discrete energy levels even though the 
potential for e = 1 is imaginary and the potential for e = 2 is upside-down. The first 
complete proof of spectral reality and positivity for H m. ^ was given by Dorey et al. in 
Ref. 0,13. 



The philosophical background of PT quantum mechanics is simply this: One of the 
axioms of quantum mechanics requires that the Hamiltonian H be Dirac Hermitian. This 
axiom is distinct from all other quantum-mechanical axioms because it is mathematical 
rather than physical in character. The other axioms of quantum mechanics are stated in 
physical terms; these other axioms require locahty, causality, stability and uniqueness of the 
vacuum state, conservation of probability, Lorentz invariance, and so on. The condition of 
Dirac Hermiticity H = W is mathematical, but the condition of PT symmetry H = H^^ = 
{PT)H{PT) (space-time reflection symmetry) is physical because P and T are elements of 
the Lorentz group. 

The spectrum of if in ([2]) is real, which poses the question of whether this Hamiltonian 
specifies a quantum-mechanical theory. That is, is the theory specified by H associated 
with a Hilbert space endowed with a positive inner product and does H specify unitary 
(norm-preserving) time evolution? The answer to these questions is yes. Positivity of the 



inner product and unitary time evolution was established in Ref. [11|, Il2| for quantum- 
mechanical systems havirig an unbroken PT symmetry (an analogous result was obtained 
by Mostafazadeh in Ref. 5^]) and in Ref. [l3| for quantum field theory. 



To demonstrate that the theory specified by the if in ([2]) is a quantum-mechanical theory, 
we construct a linear operator C that satisfies the three simultaneous algebraic equations 



ll| : = 1, \C,PT] = 0, and [C,H] = 0. Using C, which in quantum field theory is a 
Lorentz scalar [s}], we can then construct the appropriate inner product for a PT-symmetric 
Hamiltonian: {a\b) = a^^'^ ■ b. This inner product, which uses the CPT adjoint, has 
a strictly positive norm: {a\a) > 0. Because H commutes with both PT and C, H is 
self-adjoint with respect to CPT conjugation. Also, the time-evolution operator q-^^^I^^ ig 
unitary with respect to CPT conjugation. Note that the Hilbert space and the CPT inner 
product are dynamically determined by the Hamiltonian itself. 

We have explained why a PT-symmetric Hamiltonian gives rise to a unitary theory, but 
in doing so we raise the question of whether PT-symmetric Hamiltonians are useful. The 
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answer to this question is simply that PT-symmetric Hamiltonians have already been useful 
in many areas of physics. For example, in 1959 Wu showed that the ground state of a Bose 
system of hard spheres is described by a non-Hermitian Hamiltonian ^]. Wu found that 
the ground-state energy of this system is real and he conjectured that all of the energy levels 
were real. Hollowood showed that the non-Hermitian Hamiltonian for a complex Toda lattice 
has real energy levels 4l|]. Cubic non-Hermitian Hamiltonians of the form H = p^ + ix^ (and 



also cubic quantum field theories ha yirig an irnaginary self-coupling term) arise in studies 
of the Lee- Yang edge singularity 12 
[2^ 



25l . l3ll . |61| and in various Reggeon field-theory models 



In all of these cases a non-Hermitian Hamiltonian having a real spectrum 
appeared mysterious at the time, but now the explanation is clear: In every case the non- 
Hermitian Hamiltonian is PT symmetric. Hamiltonians having PT symmetry have also 
been used to describe magnetohydrodynamic systems [13, UTi] and to study nondissipative 
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time-dependent systems interacting with electromagnetic fields 

An important application of PT quantum mechanics is in the revitalisation of theories 
that have been thought to be dead because they appear to have ghosts. Ghosts are states 
having negative norm. We have explained above that in order to construct the quantum- 
mechanical theory defined by a PT-symmetric Hamiltonian, we must construct the appro- 
priate adjoint from the C operator. Having constructed the CPT adjoint, one may find that 
the so-called ghost state is actually not a ghost at all because when its norm is calculated 
using the appropriate definition of the adjoint, the norm turns out to be positive. This is 
what happens in the case of the Lee model. 

The Lee model was proposed in 1954 as a quantum field theory in which mass, wave- 
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function, and charge reorganisation could be performed exactly and in closed form 
However, in 1955 Kallen and Pauli showed that when the renormalised coupling constant is 
larger than a critical value, the Hamiltonian becomes non-Hermitian (in the Dirac sense) and 
a ghost state appears 4J] . The importance of the work of Kallen and Pauli was emphasised 



by Salam in his review of their paper [57[ and the appearance of the ghost was assumed 
to be a fundamental defect of the Lee model. However, in 2005 it was shown that the 
non-Hermitian Lee-model Hamiltonian is PT symmetric and when the norms of the states 
of this model are determined using the C operator, which can be calculated exactly and in 
closed form, the ghost state is seen to be an ordinary physical state having positive norm . 
Thus, the following assertion by Barton is not correct: "A non-Hermitian Hamiltonian is 
unacceptable partly because it may lead to complex energy eigenvalues, but chiefly because 
it implies a non-unitary S matrix, which fails to conserve probability and makes a hash of 
the physical interpretation." 

Another example of a quantum model that was thought to have ghost states, but in fact 



does not, is the Pais-Uhlenbeck oscillator model 16, 17, 18 . This model has a fourth-order 



field equation, and for the past several decades it was thought (incorrectly) that all such 



higher-order field equations lead inevitably to ghosts. It is shown in Ref. [16|] that when the 
Pais-Uhlenbeck model is quantised using the methods of PT quantum mechanics, it does 
not have any ghost states at all. 

There are many potential applications for PT quantum mechanics in areas such as par- 
ticle physics, cosmology, gravitation, quantum field theory, and solid-state physics. These 
applications are discussed in detail in the recent review article in Ref. Furthermore, 
there are now indications that theories described by PT-symmetric Hamiltonians can be 



observed in table-top experiments [48|, [54 



Having shown the validity and potential usefulness of PT quantum mechanics, one may 
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ask why PT quantum mechanics works. The reason is that CP is a positive operator, and 
thus it can be written as the exponential of another operator Q: CP = e^. The square root 
of e*^ can then be used to construct a new Hamiltonian H via a similarity transformation 
on the PT-symmetric Hamiltonian H: H = e~*^/^ife'^/^. The new Hamiltonian H has the 
same energy eigenvalues as the original Hamiltonian H because a similarity transformation 



is isospectral. Moreover, H is Dime Hermitian [51| ; PT quantum mechanics works because 
there is an isospectral equivalence between a non-Hermitian PT-symmetric Hamiltonian and 
a conventional Dirac Hermitian Hamiltonian. 

There are a number of elementary examples of this equivalence, but a nontrivial illustra- 
tion is provided by the Hamiltonian i7 in ([2]) at e = 2. This Hamiltonian is not Hermitian 
because boundary conditions that violate the norm must be imposed in Stokes wedges in 
the complex plane in order to obtain a real, positive, discrete spectrum. The exact equivalent 



Hermitian Hamiltonian is H = p'^ + Ax^ — 2Hx, where H is Planck's constant [10|, l23|, |43j . The 
term proportional to H vanishes in the classical limit and is thus an example of a quantum 
anomaly. 

Since PT symmetry is equivalent by means of a similarity transformation to conventional 
Dirac Hermiticity, one may wonder whether PT quantum mechanics is actually fundamen- 
tally different from ordinary quantum mechanics. The answer is yes and, at least in prin- 
ciple, there is an experimentally observable difference between PT-symmetric and ordinary 
Dirac Hermitian Hamiltonians. The quantum brachistochrone provides a setting for exam- 
ining this difference and provides a way to discriminate between the class of PT-symmetric 
Hamiltonians and the class of Dirac Hermitian Hamiltonians. 



III. COMPLEX CLASSICAL MOTION 

It is implicitly assumed in the derivation of the classical brachistochrone that the path of 
shortest time of descent is real. However, it is interesting that if one allows for the possibility 
of complex paths of motion, one can achieve an even shorter time of flight. In this section 
we consider a simple classical-mechanical system. Our purpose is to explain heuristically 
how extending a dynamical system into the complex domain can result in faster-than-real 
time evolution. 

To demonstrate that a shorter time of flight can be achieved by means of complex paths, 
let us consider the classical harmonic oscillator, whose Hamiltonian is 

H = p'^ + x\ (3) 

If a particle has energy E = 1, then the classical turning points of the motion of the particle 
are located at x = ±1. The particle undergoes simple harmonic motion in which it oscillates 
sinusoidally between these two turning points. This periodic motion is indicated in Fig. [1] 
by a solid line connecting the turning points. However, in addition to this oscillatory motion 
on the real-x axis, there are an infinite number of other trajectories that a particle of energy 
E can have 0]. These classical trajectories, which are also shown in Fig. [T], are all ellipses 
whose foci are located at precisely the positions of the turning points. All of the classical 
orbits are periodic, and all orbits have the same period T = 2tt. Thus, a classical particle 
travels faster along more distant ellipses. 

Now suppose that a classical particle of energy E = 1 travels along the real-x axis from 
some point x = —a to x = a, where a > 1. If the potential V{x) is everywhere zero along its 
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FIG. 1: Classical trajectories in the complex-x plane for the harmonic oscillator whose Hamiltonian 
is H = + x^. These trajectories represent the possible paths of a particle whose energy \s E = \. 
The trajectories are nested ellipses with foci located at the turning points at a; = ±1. The real line 
segment (degenerate ellipse) connecting the turning points is the usual periodic classical solution 
to the harmonic oscillator. All closed paths have the same period 27r. 

path, then it will travel at a constant velocity. However, suppose that the particle suddenly 
finds itself in the parabolic potential V{x) — when it reaches the turning point x — —1 
and that it suddenly escapes the influence of this potential aX, x — 1. Then, the time of flight 
from X = —a to x = a will be changed because the particle does not travel at a constant 
velocity between the turning points. Next, let us imagine that the potential V{x) = is 
suddenly turned on before the particle reaches the turning point at a; = — 1. In this case, the 
particle will follow one of the elliptical paths in the complex plane around to the positive 
real axis. Just as the particle reaches the positive real axis the potential is turned off, so 
the particle proceeds onward along the real axis until it reaches x = a. This trip will take 
less time because the particle travels faster along the ellipse in the complex plane. 

We have arrived at the surprising conclusion that if the classical particle enters the 
parabolic potential V{x) = x^ immediately after it begins its voyage up the real axis, its 
time of flight will be exactly half a period, or vr. Indeed, by travelling in the complex plane, 
a particle of energy E — \ can go from the point x = —a to the point a; = a in time 
TT, no matter how large a is. Evidently, if a particle is allowed to follow complex classical 
trajectories, then it is possible to make a drastic reduction in its time of flight between two 
given real points. 
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IV. HERMITIAN QUANTUM BRACHISTOCHRONE 

The quantum brachistochrone problem, as described briefly in Sec. HI is similar to the 
classical counterpart except that the optimisation takes place in a Hilbert space. Specifically, 
we are given a pair of quantum states, an initial state and the final state \iPf)i and 
we would like to find the one-parameter family of unitary operators {Ut} that achieves 
the transformation \'4^f) = Ut\ipi) in the smallest possible time t. Since a one- 

parameter family of unitary operators can be formed in terms of a Hermitian operator H 
as Ut = exp{—iHt/fi), the problem is equivalent to finding the Hermitian operator H that 
realises the transformation {"ipp) in the shortest possible time. 

The Hermitian operator H can be thought of as representing the Hamiltonian, so the 
quantum brachistochrone problem is equivalent to finding the optimal Hamiltonian H satis- 
fying exp{—iHt / H)\ip i) = I'ipF)- However, it is intuitively clear that if we are allowed to have 
access to an unbounded energy resource, then the time required for the relevant transforma- 
tion, irrespective of whether the Hamiltonian is optimal or not, can be made arbitrary small. 
Hence, for a quantum brachistochrone problem to possess a nontrivial solution, some form 
of constraint is needed. The simplest constraint is to assume that the energy is bounded so 
that the difference between the largest and the smallest energy eigenvalues has a fixed value 

-^max -^min (4) 

A short calculation shows that if the Hamiltonian H is bounded, then (i) the standard 
deviation of the energy is bounded according to 

AiJ < |(-Emax — -E'min), (5) 

and (ii) the state with maximum energy uncertainty is (|-Emax) + \ Emm))/ V^- It follows that 
the energy constraint (jl]) is equivalent to a constraint on energy uncertainty. 

The brachistochrone problem of this type has been analysed recently and a solution was 
obtained by means of a variational method 2^ . It has also been solved in terms of a more 



elementary approach making use of the geometry of quantum state space [2C|. We shall 
follow closely the latter approach here. 

Let us now state the simplest form of the quantum brachistochrone problem: We have 
a quantum system represented by an iV-dimensional Hilbert space Ti and a prescribed pair 
of states and {ipp) on Ti. The problem is (a) to find the Hamiltonian H satisfying the 
constraint (jlj) such that the unitary transformation exp{—iHt / fi)\iljj) = \iPf) is achieved in 
shortest possible time and (b) to find the time required to realise such an operation. 

A little geometric intuition allows us to find the solution to this problem with minimum 
effort. Recall that the time required for transporting a state along a path in Ti is given by 
the distance divided by the speed. Hence, all we have to do is first to identify the shortest 
path and measure its length and then to allow the state to evolve along the path with the 
greatest possible speed without violating the constraint (jl]). 

In quantum mechanics the notion of distance is closely linked to the notion of transition 
probability (21I . 4^ . In particular, by looking at the transition probability between neigh- 



bouring states we can derive the expression for the metric on the space of quantum states. 
This allows us to measure distances between states. The idea can be sketched as follows: 
Consider a state {ip) in Ti. and a neighbouring state {ip) + \dtp). The transition probability 
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between these states is 



cos 



2 (IaA = V\V"I \'-^'^i;i¥^/\y|V|y/ "T l^H^/J /p\ 



(^i^)((^i + (d^i)(i^) + id^)) 

where ds defines the fine element on the space of pure states. By using 



cos 



^ (Id.) - 1 - ]ds\ (7) 
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expanding the right side of ([6]), and retaining terms of quadratic order, we find that the fine 
element is 

2 _ , (^|^)(dV;|d^)-(^|dV;)(d^|^) 

im' ■ ^ ^ 

This line element is known in geometry to arise from the Fubini- Study metric and it can 
be used to measure the distance of the shortest path joining a pair of points on the space of 
pure quantum states. 

If the Hilbert space is two dimensional, then a generic normalised state vector can be 
expressed in the form 



cos \6 



sin 16* e^"^ 



(9) 



A short calculation then shows that the Fubini-Study line element ([8]) reduces in this case 
to the expression 

ds^ = d^2^sin2^d02, (10) 

which we recognise as the line element on the Bloch sphere S. (The Bloch sphere is the 
state space of two-level systems.) 

In the case of an n- dimensional Hilbert space 7i, if we are given a pair of distinct states 
and \iPf)i then the linear span of these two states forms a two-dimensional subspace of 
Ti. It should be intuitively clear that the shortest path joining I?/'/) and \'4'f) should lie on this 
two-dimensional subspace. Thus, irrespective of the dimensionality of Ti the solution to our 
quantum brachistochrone problem can be obtained by analysing the two-dimensional subspace 
spanned by {ipj) and \iPf)- Even when we restrict our attention to this subspace, there still 
remain infinitely many unitary orbits that realise the transformation {"ipF) = Ut\ipi). 

However, since the two-dimensional state space is just the Bloch sphere S endowed with the 
spherical metric (fTOj) . we see that there is a unique great circle arc that joins {ipi) and \iPf) 
on S. (This assumes, of course, that {ipi) and \iPf) are not antipodal points of the sphere. 
Otherwise, there are infinitely many such paths.) In this way we have identified the shortest 
path joining and I^Pf)- The shortest distance Smin between these two points of S is thus 
given by 

..in = 2arccosf^MfcLV (11) 
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This result can also be obtained by integrating the line element (fTOl) along the great-circle 
arc on S. 

Having obtained the distance of the shortest path we proceed to find the maximum speed 
at which the state can evolve unitarily. For the evolution of the state we must consider the 
general Schrodinger equation, but we also need to express the equation in the correct form. 
This is the so-called modified Schrodinger equation 

In this equation the mean- adjusted Hamiltonian H is given by 

H = H-{H), (13) 

where 

Note that {H) = and that according to fll2p the tangent vector ^\psit) is everywhere 
orthogonal to the direction of the state [45]. Since the energy expectation (H) depends on 
the state the modified Schrodinger equation appears to be nonlinear. However, this is 
not the case. The point is that the expectation value of the Hamiltonian is a constant of 
the motion under the Schrodinger dynamics. Thus, given the initial state we calculate 
(H) in this state and subtract this number from the Hamiltonian. Since the Hamiltonian in 
quantum mechanics is defined only up to an additive constant, this modification does not 
alter the physics in any way. It is worthwhile noting that the modified Schrodinger equation 
ffT^ is canonical and reduces to the standard eigenvalue problem when the state lipt) is time 
independent without one having to evoke the correspondence principle ^2l\. 

If the initial state vector lipj) is normalised, then the evolution f[T^ preserves the norm. 
It follows that {ipldip) = 0. Since the speed v of quantum evolution is given hj v = ds/dt, 
we find from and that 

v' = ^, mn - (iJ))V*) = ^ {HiH - {H)f\i^i)- (15) 

This shows that the speed of quantum evolution is given by the energy uncertainty. The 
expression (fTSl) for the speed of quantum evolution is known as the Anandan-Aharonov 
relation fl|]. Since we know from (JSj) that under the constraint (Hj) the energy uncertainty 
A/7 is bounded by |ci;, we find that the maximum speed of quantum evolution is given by 

?^max = ^- (16) 

By using the results in (fTTl) and (fT6|) we deduce that the minimum time required for realising 
the unitary transportation \ipi) l^pp) = Ut\ipi) is given by the ratio Smin/'^max- In particu- 
lar, if and {ipp) are orthogonal, then they correspond to antipodal points on the Bloch 
sphere S, and we have Smm = tt. In this case, the minimum time required to orthogonalize 
the state (that is, for the state to evolve into a new state that is orthogonal to the original 



state) is known as the passage time rp [19|, l58|. The passage time is explicitly 
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FIG. 2: Optimal Hamiltonian for quantum state transportation. The two-dimensional complex 
Hilbert space spanned by the initial state IV'/) and the final state \^f) can be visualised in real 
terms as a Bloch sphere S. The two states and \^f) can then be identified as a pair of points on 
S. Assuming that these two points are not antipodal, there exists a unique great circle arc joining 
these two points, which determines the shortest path joining the two states. The optimal way of 
unitarily transporting \ipj) into \iPf) is therefore to rotate the sphere along the axis orthogonal to 
the great circle. The axis of rotation then specifies two quantum states: \E-) and These 
states are the eigenstates of the optimal Hamiltonian H. 

The passage time f|T7|) provides the bound in Hermitian quantum mechanics for transporting 
a state into an orthogonal state, and is sometimes referred to as the Fleming bound ^] . 

The ratio Smm/'^max gives the solution to part (b) of our quantum brachistochrone prob- 
lem. To solve part (a), that is, to find the optimal Hamiltonian, we argue as follows: Since 
the problem is confined to a two-dimensional subspace of Ti, the solution can be obtained by 
elementary trigonometry on the Bloch sphere S. The key idea is to recall that the shortest 
path joining and \ipp) is a great circle arc on S. Without loss of generality, we can 
assume that {ipj) and {ipp) lie on the equator of S with respect to a suitable choice of axis. 
Then, the unitary motion along the shortest path can be generated by the rotation of S 
along this axis. Since the eigenstates of the Hamiltonian H that generates such a rotation 
correspond to the antipodal points along this axis, the pair of states \ipi) and \iPf) can both 
be expressed as equal superpositions of the eigenstates of H with the relative phase shifted 
by Sjnin- Writing \E+) and \E_) for the normalised eigenstates of H and using a = Smin/2, 
we can express the initial and the final state in the form 

^(|E_)+e-|E+)) = |^,) and j-^{\E_) + e'-\E^)) = \^f). (18) 

Solving these equations for \E^) and \E_), we obtain 
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and 

These are the eigenstates of the optimal Hamiltonian H that generates the unitary motion 
iV'-f) ~^ I'^f) = Utlipi) along the shortest path. The eigenvalues of the optimal H can be 
arbitrary as long as the condition (jlj) is satisfied. Without loss of generality, we may assume 
H be trace free, and we obtain the solution to part (a) of the quantum brachistochrone 
problem: 

H =\uj\E+){E+\-\uj\E^){E4. (21) 

This is the 'minimal' solution to the problem in the sense that H acts only on the two- 
dimensional subspace of Ti while leaving the rest of Ti unperturbed. 

As a special example, consider the problem of a spin-flip, that is, turning a spin-up state 
into a spin-down state unitarily. In this case the initial and the final states can be written 

as 

\i^i)=(]\ and \iJF)={^\ (22) 



.0/ \l, 

in the spin-z basis. Substituting these into (IT^ and (12(1 . we find that the eigenstates of the 
optimal Hamiltonian are 

because in this case we have a = 7r/2. Substituting this result into ( 12T1) yields 

for the optimal Hamiltonian. By using the relation 

where n is a unit vector and 

^i=(lo)' ^^=(z7)' ^^=(o-l) (26) 

are Pauli matrices, we obtain the expression for the optimal unitary operator: 

/ cos(i) -isin(i)\ 
Ut={ . (27) 

V-isin(i) cos(i); 

It follows at once that the optimal unitary orbit {ipt) = Ut\ipi) is given by 

\^^)=i '"^fi 1- (28) 
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We find that the first time at which \'il)t) reaches IV'f) is given by the condition ujt/2fi = 7i/2, 
that is, when t = rp, where rp is the passage time given in (fT71) . 

We have seen how the simplest form of a quantum brachistochrone problem can be 
solved in Hermitian quantum mechanics by considering a two-dimensional Hilbert subspace 
combined with elementary geometric constructions on it. In a more general situation the 
unitary motion may be constrained further so that the optimal Hamiltonian (|2T]) may not be 
implementable. For example, the constraint may enforce the path of the unitary evolution 
to lie in a three- rather than in a two-dimensional subspace. To determine what happens 
let us work out the passage time for this example. 

Since in this case the minimal solution H to the brachistochrone problem is a three- 
dimensional matrix, we can express the initial state {ipi) in terms of the three eigenstates of 
H according to 

lijji) = cosa\Ei) + sin a cos /3e^'^\Ej) + sina sin/3e"^|i?fc), (29) 

where a, [3 are angular coordinates, 0, ip are phase variables, and we assume that Ei < 
Ej < Ek- If a unitary operator Ut transforms this state into an orthogonal state, then the 
condition 

cos^ a + sin^ a cos^ /^e'^'^J"^/^ + sin^ a sin^ fj^-i^MT/n ^ q ^gg^ 

must be satisfied, where ujji = Ej — E^ and uj^i = E^ — E.-^. To render the analysis more 
tractable, we simplify this constraint by assuming that a = (3 = tt / A. Then flSU]) implies 
that a necessary condition for the state to evolve into an orthogonal state is given by 
the relation 

ujki 2m — 1 
UJji 2n — l' 

where m, n are natural numbers such that m ^ n. 

Thus, the solution to the brachistochrone problem must be such that the eigenvalues of 
H satisfy condition (1311) as well as the constraint -Emax ~ -^'min ^ ^- Assuming that these 
constraints are indeed satisfied, the initial state evolves into an orthogonal state IV'f)- The 
first time that \ipi) evolves into a state orthogonal to lipi), in particular, is given by 

UJji UJki 

However, since in this case Ut\ilJi) does not describe a geodesic path, T will be larger than 
Fleming's passage time rp given in f|T71) . Indeed, without loss of generality we may set 
Ei = 0. Then, it is straightforward to verify that T = y/Erp. This follows from the fact that 
under the constraint Uki = Sujji that comes from ( l32|l . the squared energy dispersion in the 
state (EHl) with a = (3 = n/4 is given by AH"^ = 

V. NON-HERMITIAN QUANTUM BRACHISTOCHRONE 

We have seen how the solution to the simple brachistochrone problem can be obtained 
in the Hermitian quantum theory. What happens if we extend the quantum theory into 
the complex domain by looking at a PT-symmetric theory? We saw earlier that in classical 



(31) 



(32) 
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mechanics if we were to allow for a complex path interpolating a pair of real points of the 
coordinate space, then it is possible (at least mathematically) to transport a particle across 
a large distance in virtually no time. It turns out that an analogous situation emerges in 
the PT-symmetric theory. Here we present a simple algebraic calculation of the optimal 
evolution time from an initial state to a final state by using a simple 2x2 Hamiltonian. 
As we have remarked above, the 2x2 model suffices to cover general cases because in the 
case of our simple brachistochrone problem the solution is found on the two-dimensional 
subspace of the Hilbert space spanned by the initial state and the final state \iPf)- In 
the case of a PT-symmetric Hamiltonian the variational approach gives a more direct way 
to handle the brachistochrone problem. Thus, we shall first briefiy revisit the Hermitian 
case but expressed in the variational formalism and then we will compare the result to its 
PT-symmetric counterpart. 

A. Hermitian 2x2 matrices 

We choose a basis so that the initial and final state vectors take the form 

I^/)=(J) and I^f) = (^), (33) 

where the condition that \iPf) be normalised is |ap + |6p = 1. The most general 2x2 
Hermitian Hamiltonian is 

(^Ji. '^''J'') (r, u, ^ real). (34) 

For this Hamiltonian the eigenvalue constraint (jl]) takes the form 

= {s-uf + Ar\ (35) 

To find the optimal Hamiltonian satisfying this constraint, we rewrite H as a linear 
combination of Pauli matrices: 

H = l{s + u)t + luja-n, (36) 

where ^ 

n = — (2r cos 6', 2r sin 6', s — m) (37) 

is a unit vector. Then by use of the identity ( l25l) the evolution equation \iPf) = e~^^'^^^\'ijji) 
can be expressed in the form 



-^i{s+u)t/h 



COS g - i^ sin g 

2n u) 2n 



(38) 



The second component of this equation gives |6| = ^sin||, which allows us to find the 



required time of evolution: 



2^ . ^I&l 

t = — arcsm — — . (39) 
UJ 2r 
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We must now minimise the time t over all r > while maintaining the energy constraint 
in (135|) . This constraint tells us that the maximum value of r is ^u. At this value we have 
s = u. Because H can be made trace free, we can set s = m = 0. The variable 9 in fl36l) does 



not affect the eigenvalues, so we may set 9 = tt. Then we recover the optimal Hamiltonian 
obtained in fl2^ . As regards the minimum evolution time r we have 

^ ^ 2^arcsin|&| ^ ^^^^ 

UJ 

In the special case for which a = and 6 = 1 so that and \iPf) are orthogonal, we 
recover the passage time r = rp = ttH/uj, the smallest time required for a spin flip. 

Although the form of the result in fHUl) resembles the statement of the uncertainty prin- 
ciple, it is merely the statement indicated above that rate xtime= distance; the maximum 
speed of evolution is given by AH, and the distance between \ipi) and lipp), assuming they 
are normahsed, is given by 2 arccos(| (t/'fIV'/) I)- Since {{ippl'ipi)] = \a\ and \a\ = -y/l — 
we obtain (HOj) from the relation 

arccosA/l — = arcsin|6|. (41) 



B. Non-Hermitian 2x2 matrices 

We now show by direct calculation that for a PT-symmetric Hamiltonian, r can be 
arbitrarily small. This is because a PT-symmetric Hamiltonian whose eigenvalues are all real 
is equivalent to a Hermitian Hamiltonian via H = e~'^^'^He^^'^, where Q is Dirac Hermitian. 
The states in a PT-symmetric theory are mapped by e"*^/^ to the corresponding states 
in the Dirac Hermitian theory. But, the overlap distance between two states does not 
remain constant under a similarity transformation. We can exploit this property of the 
similarity transformation to overcome the Hermitian lower limit on the time r. (The detailed 
calculation is explained in Ref. [14].) 

We consider the general class of PT-symmetric 2x2 Hamiltonians having the form 

^=('f^,-i.) (r, ^real). (42) 

The time reversal operator T performs complex conjugation and the parity operator in this 
case is given by 



1 

1 



(43) 



The two eigenvalues 



E± = rcos9 ± a/s^ - sin^ 9 (44) 

are real if > sin^ 9. This inequality defines the region of unbroken PT symmetry. The 
unnormalised eigenstates of H are 

ia/2 \ f ip-ia/2 



= ( :-W2 and \EJ) = _ ) , (45) 
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where a is given by sin a = (r/s) sin^. Note that the condition of unbroken PT symmetry 
of H in (H2|) imphes that a is reaL The C operator required for defining the Hilbert space 
inner product is 

C=^P^\^« V (46) 

cosa V 1 -ismay ' 

It is easy to verify that the CPT norms of both eigenstates have the value \/2 cos a. 
To calculate r we express the Hamiltonian H in fH2|) as 



= (rcos^)l + iu;a-n, (47) 

where 

2 

n = — (s,0,irsin6') (48) 

UJ 

is a unit vector. The energy constraint requires that the squared difference between energy 
eigenvalues is 

The positivity of tu^ is ensured by the condition of unbroken PT symmetry. Notice that fH^ 
differs from (!35l) by a sign. We can think of (149!) as being hyperbolic in character, while (135|) 
is elliptic in character. The technical advantage of the constraint in (149!) is that because of 
the minus sign both terms on the right side can become large without violating the condition 
that UJ be fixed. We will see that it is this fact that allows the non-Hermitian Hamiltonian 
H in (H2l) to achieve faster-than-Hermitian time evolution. 

To determine r we write down the PT-symmetric time-evolution equation in vector form: 

/ -I \ ^-itrcose/h ( cos(|| - a) \ 
cos a 

In particular, consider the pair of vectors used in the Hermitian spin-flip case as in (l22l) . 
Observe that the evolution time t needed to reach \il)p) = (°) from = (g) is given by 

t = (^^HiO^. (51) 

UJ 

Optimising this result over allowable values for a as a approaches ^tt, the optimal time r 
tends to zero, a dramatic change from the Hermitian result in (1T71) ! Note, however, that the 
two vectors in (122!) are not orthogonal with respect to the CPT inner product. This is the 
reason that the Fleming bound in (1T7|) is not violated. 



VI. EXTENSION OF NON-HERMITIAN HAMILTONIANS TO 
HIGHER-DIMENSIONAL HERMITIAN HAMILTONIANS 

We have seen how a quantum state can be transported unitarily into another state in 
arbitrary short time by using a bounded Hamiltonian if we allow for a complex path in- 
terpolating them in the space of unitary motions. Can such an operation be implemented 
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in practice? If the answer is affirmative, then the imphcation is immense in quantum in- 
formation, computation, cryptography, and other related fields. For example, if a quantum 
computer were to exist, then solutions to difficult optimisation problems can in principle be 
found in arbitrary short time, and this in turn would have important implications in society 
as a whole. 

A gedanken experiment was proposed in Ref. [3] to realise this effect in a laboratory. 
The setup is as follows: we use a Stern-Gerlach filter to create a beam of spin-up electrons. 
The beam then passes through a black box containing a device governed by a PT-symmetric 
Hamiltonian that flips the spins unitarily in a very short time. The outgoing beam then 
enters a second Stern-Gerlach device that verifies that the electrons are now in spin-down 
states. In effect, the black box device is applying a complex magnetic field B: 

B = {s,0,irsme). (52) 

If the field strength has sufficiently large amplitude, then spins can be flipped in virtually no 
time because the complex path joining these two states is arbitrarily short without violating 
the energy constraint in fH9|) . We emphasise that the field strength can be made large 
without violating the energy constraint (jlj) is a consequence of the hyperbolic representation 
in (Uni). 

The arbitrarily short alternative complex pathway from an up state to a down state, 
as illustrated by this thought experiment, is reminiscent of the short alternative distance 
between two widely separated space-time points as measured through a wormhole in general 
relativity. This comparison is of course controversial, and it has subsequently motivated 
much research and it has generated a lively debate in the literature [l, [H, 33, 35, [1^, H^, 
52, 53, 5^. We emphasise that the entire package of flipping the spin is not realised by 
a unitary operation. This follows from the fact that PT-symmetric quantum theory is 
unitary, and as such it respects the Fleming bound (1171) applicable to all unitary theories 



14l |. The point is that the 'black box' scheme described above actually consists of three 
regimes: (i) the preparation of a spin-up state in the Hermitian setup, (ii) the fast unitary 
motion to ffip the spin using a PT-symmetric Hamiltonian, and (iii) the recovery of a spin- 
down state in the Hermitian setup. Thus, the operation is locally unitary, but the switching 
between Hermitian and PT-symmetric description is not unitary. This three-step process 
of switching Hamiltonians is analogous to the classical procedure described in Sec. IIIII for 
obtaining faster-than-real time evolution. Recall that in the classical case we were able to 
transport a finite-energy particle across a large distance in a short time by switching the 
potential through which it was travelling. Note that in the classical case there is no question 
of violating unitarity because the particle does not get lost. 

The question of unita rity and faster-than-Hermitian time evolution has been reexamined 
in more detail in Ref. 35, 53] by means of a geometric approach and also in Ref. 0], 
where a more general class of non-Hermitian Hamiltonians that are not necessarily PT 
symmetric are considered. In particular, Mostafazadeh has emphasised the role of quantum 
observables in such an experiment; the spin operator in Hermitian quantum mechanics 
cannot be interpreted as a spin operator in the PT-symmetric counterpart, thus leading 
to ambiguities regarding the physical interpretation of the gedanken experiment described 
above. 

An intriguing alternative proposal for an implementation of the fast spin flip has been 
made more recently by Giinther and Samsonov ^] . The idea is to embed the problem into 
a Hermitian setup represented by a higher-dimensional Hilbert space. Specifically, take our 
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PT-symmetric Hamiltonian H in (H2|) . The eigenstates of H are not orthogonal with respect 
to the Hermitian inner product. Since H is not Hermitian, its Hermitian conjugate defines 
a new matrix H'^ . The eigenstates of W thus also define a pair of nonorthogonal states in 
the Hermitian theory. When these four states are suitably normalised, they can be used 
to form an over-complete basis set in the Hermitian two-dimensional Hilbert space. Such 
an over-complete set of basis is also known as a positive operator- valued measure (POVM), 
commonly used in the analysis of quantum information theory. A key idea is that such a 
basis can be embedded in a higher- dimensional Hilbert space to form an orthogonal basis by 
using the Naimark dilation [40] . A Hermitian Hamiltonian can then be constructed — in this 
case a 4 X 4 matrix — such that its eigenstates are precisely the four states thus obtained. 
Using this Hamiltonian it is possible to construct a standard unitary motion in such a way 
that the induced motion obtained by the projection onto the two-dimensional subspace is 
characterised by the PT-symmetric motion ( l50l) . 

In this way, Giinther and Samsonov were able to show that the fast spin fiip can in 
principle be realised in the standard Hermitian quantum mechanics by a combination of a 
unitary motion and a projection in a larger-dimensional Hilbert space. In practical terms 
this means that one should couple the spin to an auxiliary particle (this can be done either 
by a projection or by a unitary operation), apply a unitary evolution in the larger Hilbert 
space of the combined system, and finally project out the auxiliary particle to recover the 
spin in the transported state. The net effect of such an operation can then be characterised 
by (!50l) . The Fleming bound is apparently violated due to the general fact that when a 
unitary motion is projected to a subspace of a Hilbert space, the resulting dynamics need 
not respect laws of unitarity. It would be of considerable interest to find out whether the 
Giinther and Samsonov scheme can actually be implemented in a laboratory, and if not, 
what might be the difficulty preventing the violation of the Fleming bound. 
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